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Abstract 

Combining and expanding on work from previous publications, a model for the evolution of 
ultrarelativistic heavy-ion collisions at the CERN SPS for 158 AGeV beam energy is presented. 
Based on the assumption of thermalization and a parametrization of the space-time expansion of 
the produced matter, this model is able to describe a large set of observables including hadronic 
momentum spectra, correlations and abundancies, the emission of real photons, dilepton radiation 
and the suppression pattern of charmonia. Each of these obervables provides unique capabilities to 
study the reaction dynamics and taken together they form a strong and consistent picture of the 
evolving system. Based on the emission of hard photons, we argue that a strongly interacting, hot 
and dense system with temperatures above 250 MeV has to be created early in the reaction. Such 
a system is bound to be different from hadronic matter and likely to be a quark-gluon plasma, and 
we find that this assumption is in line with the subsequent evolution of the system that is reflected 
in other observables. 



1 Introduction 

Lattice simulations (e.g. [1]) predict that QCD undergoes a phase transformation at a temperature 
Tc ~ 150 — 170 MeV [2,3] from a confined hadronic phase to a phase, the quark-gluon plasma (QGP), in 
which quarks and gluons constitute the relevant degrees of freedom and the chiral condensate vanishes. 
Experimentally, this prediction can only be tested in ultrarelativistic heavy-ion collisions. However, 
finding evidence, and ultimately proof, for the creation of a quark-gluon plasma faces several difficulties. 
Arguably the greatest challenge is to link experimental observables to quantities measured on the lattice. 

For a large set of observables, the evolution of the expanding medium is a key ingredient for their 
theoretical description. It is also the place where results from lattice QCD fit in: assuming a thermalized 
system is created, its evolution is governed by the equation of state (EoS). Hence, in order to test 
the lattice QCD predictions, one has to start with this assumption and show that it leads to a good 
description of the experimental data. 

It is the purpose of the present paper to show that such a description can be achieved in the case of 158 
AGeV Pb-Pb and Pb-Au coUisions at the CERN SPS. We summarize and expand on the results of [4-10] 
by developing a framework based on the EoS determined in lattice QCD. Using eikonal calculations 
to constrain the initial state and measured bulk hadronic properties to establish the final state, we 
parametrize the evolution in between in a way that is motivated by hydrodynamical calculations. We 
then use this parametrization to calculate other observables not connected to the bulk hadronic matter 
such as dilepton and photon emission and charmonium suppression. In all cases, we find good agreement 
with the data with the same set of model parameters. We explicitly discuss how each observable reflects 
underlying scales of the evolution, to what degree information from a particular observable constrains 
alternative evolution scenarios and what we can learn about the relevant degrees of freedom in the 
medium. 



2 The fireball model framework 



2.1 Expansion and flow 

We do not aim at a microscopical description of the expansion but instead at an effective parametrization 
based on available data. In order to simplify computations, we introduce several assumptions which lead 
to an idealized picture of the expansion process. The aim of the model is to demonstrate that there is one 
set of scales characterizing the expansion that is in reasonable agreement with a large set of observables, 
not to provide a detailed description of every observable. 

Our fundamental assumption is that an equilibrated system is formed a short proper time tq after the 
onset of the collision, which subsequently expands isentropically. As soon as the mean free path of 
particles in the medium exceeds its dimensions, kinetic decoupling occurs at a proper time r/. The 
matter then ceases to be in thermal equilibrium and we assume that no significant interactions occur 
later. 

During the expansion phase, we posit that the entropy density of the system can be described by 

s{T,7^s,r) = NR{r,T) -HiT^^.T) (1) 

with T the co-moving proper time of a given volume clement and 77^ = ^ ^^^(f^f) the spacctimc rapidity. 
i?(r, r) and H[r]s,T) are two functions describing the evolving shape of the distribution, and is a 
normalization factor chosen such that the total entropy is = / d^rsij, r]s,r). In the above expressions 
we have neglected angular asymmetries for collisions with a finite impact parameter, we expect their 
influence to be small as long as we do not discuss very peripheral collisions or observables, such as elliptic 
flow, with an explicit angular dependence. 

In oder to simplify our calculations, we choose R{r, r) and H{r]s,T) as box profiles R{r, r) = 0{R{t) — r), 

H{ris,T) = 9{Hs(t) — rjs)0{ris) + 0{ris — Hs{T))6{~ris). For this choice, thermodynamical parameters 
become functions of proper time r only. We will later argue that this simplified choice of the distribution 
function gives a fair representation of the physics at midrapidity, where most observables are measured, 
but fails to reproduce the conditions at forward rapidities. For bulk observables, like particle numbers, 
obtained by integrating over the whole volume, the detailed choice is irrelevant or reduced to a second 
order effect due to detector acceptance. 

Thus, the expansion is governed by the scale parameters R{t) and H{t). Their growth with r is a 
consequence of collective flow of the thermalized matter. For transverse flow we assume a linear relation 
between radius r and transverse rapidity p = atanh v_\_{t) = r/Rcir) ■ Pcir) with Pc{t) = atanh a±T. 
The longitudinal flow profile is dictated by the requirement that an initially homogeneous distribution 
of matter remains homogeneous at all times. 

We start with the experimentally measured width of the rapidity interval of observed hadrons 2?]^'^°"' 
at breakup. From this, we compute the longitudinal velocity of the fireball front at kinetic freeze-out 
^front -yy-g j^q^ require the initial expansion velocity Uq°"* to coincide with u^™"* but instead allow 
for a longitudinally accelerated expansion. This implies that during the evolution rj = rjs where rj is the 
momentum rapidity, is not valid . Initially, we characterize all possible trajectories of volume elements 
by a parameter c such that Vq{c) = CTq""*, c e [—1, 1]; hence we start with an expansion that is boost- 
invariant over the rapidity interval —atanh Vq°'^^ < ?7 < atanh Vq°'^^. For the final state, we require 
atanh w^™"* = t]^™"*, i.e. the distribution has to fill the experimentally observed rapidity interval. 

In order to match initial and final state, an acceleration term has to be introduced. We assume that 
the acceleration is driven by the local conditions in and around a volume element (as it would be in 

hydrodynamic calculations), hence it has to be a function of proper time. To maintain a spatially 
homogeneous distribution of matter at all times, the acceleration furthermore must correspond to a 
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re-scaling of the velocity field. Therefore, if Uzir) is the acceleration acting on the fireball front, the 
trajectory characterized by c feels the acceleration cazir). This implies 



c-v^ + cCaz{T')dT' 



y/l^v^{c,r{t',z{t'Wdt' with v%c,t)^ \ / (2) 

During the evolution, the volume element moves the distance 

z{c,t)= I v\c,T{t,z{t)))dt (3) 

in the cm. frame. We solve the set of equations numerically for all values of c by integrating the 
trajectory forward in time until a fixed r in Eq. is reached. The resulting pairs (t(c),z(c)) define a 
curve of constant proper evolution time. We can compute the length of this curve using the inavriant 
volume formula V — J da^u^ with an element of the freeze-out hypersurface and the four- velocity 
of matter. To good accuracy, volume elements lie on the curve r = const. = ^J t{c)'^ — z{cY and the 
fiow pattern can be approximated by a linear relationship between rapidity 77 and spacetime rapidity rjs 
as 77(775) — 7;^'°"* 773/7;^''°"* where 77^''°"* is computed on the c = 1 trajectory. In this case, the longitudinal 
extension in the proper volume formula yields 

with 77g™"*(T) the spacetime rapidity of the accelerating cylinder front and = 7;^™"* (r) /Ty^™"* (r) the 
mismatch between spacetime rapidity and rapidity introduced by the accelerated motion. This is an 
approximate generalization of the boost-invariant relation L{t) = 277^''°"*t which can be derived for non- 
accelerated motion. The volume of the fireball at given proper time can then be found asV — ttR{t)'^L{t) 
where we have neglected relativistic effects in the transverse dynamics (we have estimated that those 
are on the order of 10% at maximum for the volume and thus well within the parametric uncertainties). 
We illustrate these concepts in Fig. ^ 




Figure 1: Left: Different scenarios for the longitudinal fireball expansion, all tuned to produce the 
same final rapidity distribution. Shown are rapidity rj and spacetime rapidity rjs for a standard boost- 
invariant scenario without long, acceleration (solid), for complete initial stopping and re-expansion (77: 
dash-dotted 7/s: dash-dot-dotted) and and intermediate scenario (77: dotted 77^: dashed). Right: 77 as a 
function of 77s at freeze-out time r/ for the SPS scenario described in this paper. 
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2.2 Parameters of the expansion 



We now address the proper time dependence of the two scales R{t) and H{t). The parameters entermg 
^frontj-^-j (^^^^ i^j-j determined using the ansatz = ■ p{t)/€{t) which allows a soft point in 
the EoS where the ratio p/e gets small to influence the acceleration pattern. Cz and Vq°'^*' in ^ are then 
model parameters governing the longitudinal expansion and fitted to data. For the radial expansion we 
use 

R{r) =Ro+CT r dr' C dr"^^ (5) 

Jro Jro ) 

The initial radius Rq is taken from overlap calculations. This leaves a parameter ct determining the 
strength of the transverse acceleration which is also fitted to data. The final parameter characterizing 
the expansion is its endpoint given by Tf, the breakup proper time of the system. It is determined by 
the condition that the temperature T{t) drops below the freeze-out temperature Tp. 

2.3 Thermodynamics 

We calculate the total entropy So by fixing the entropy per baryon from the number of produced 
particles per unit rapidity and the number of participant baryons [4, 11]. Assuming isentropic expansion 
the entropy density s at a given proper time is then determined by s = Sq/V{t). 

We describe the EoS in the partonic phase by a quasiparticle interpretation of lattice data which has been 
shown to reproduce lattice results both at vanishing baryochemical potential /is [12] and finite fiB [13] 
and is able to extrapolate the lattice results obtained with quark masses m ^ T down to physical masses. 

For a thcrmalized hot hadronic system near the phase transition, lattice computations of the EoS offer 
little guidance since the quark masses and hence the hadron masses are still unphysically large. On the 
other hand, many of the hadronic states relevant close to the phase transition are poorly known and 
a first principle calculation of the thermodynamics of an interacting hadron gas near this point faces 
serious difficulties. We circumvent the problem by calculating thermodynamic properties of a hadron 
gas at kinetic decoupling where interactions cease to be important and an ideal gas is presumably a valid 
description. Determining the EoS at this point using the statistical hadronization described in section 
121 we choose a smooth interpolation between decoupling temperature Tp and transition temperature Tc 
to the EoS obtained in the quasiparticle description. 

Using the EoS and s(t), we compute the parameters p(r),e(r) and T(r). Since the ratio p(T)/e(r) 
appears in the expansion parametrization, we solve the model self-consistently by iteration. 

2.4 Solving the model 

Several of the model parameters are fixed by measurements, such as the total entropy 5*0 or the final 
rapidity interval 77^°"*. Others are calculated, such as the initial transverse overlap radius i?o or the 
number of participant nucleons A^part which are obtained by eikonal calculations. Only the initial rapidity 
interval ?7^™"'(to) (via c^), the final transverse velocity v±{Tf) (via cp) and the freeze-out temperature 
Tp are fit parameters and adjusted to the hadronic momentum spectra and HBT correlations. Thus, 
total multiplicity is governed by the input Sq, the relative multiplicity of different particle species is 
calculated in a statistical hadronization framework (see section |3} whereas the shape of the transverse 
mass spectra and the HBT radii determine the transverse and longitudinal expansion parameters and the 
freeze-out temperature. The equilibration time tq cannot be obtained from fits to hadronic observables 
which reflect the final state; for the time being we choose ro = 1 fm/c and discuss variations of this 
parameter later. 
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2.5 Particle emission and HBT 



We calculate particle emission throughout the whole lifetime of the fireball by evaluating the Cooper- Frye 
formula 



9 



da^p^ exp 



T, 



f 



— d'*^xS{x,p) 



(6) 



with the momentum of the emitted particle and g its degeneracy factor from the fireball surface. For 
emission hypersurfaces with a spacelikc normal we include a 0-function preventing emission from the 
surface back into the system. Note that the factor da^^p^ contains the spacetime rapidity rjs and the 
factor p^u^^ the rapidity rj. Since these are in general not the same in our model, the analytic expressions 
valid for a boost-invariant scenario [15] do not apply. 

Pionic HBT correlation radii are obtained using the common Cartesian parametrization for spin par- 
ticles 

C{q, K)-l = exp [-qlRl,,{K) - qlRl^M) - 9? <ng(^) - 2(7o9ii?^i(if )] (7) 

(see e.g. [16, 17] for an overview and further references) for the correlator. Here, K ^ ^(pi +P2) is the 
averaged momentum of the correlated pair with individual momenta pi,p2 and q = (pi — P2) the relative 
momentum. The transverse correlation radii i?out,sidc, long follow from the emission function as 



{y 



{{x - pjf) 



with P± the transverse velocity of the emitted pair, x^ 



7->2 

-'''long 

„,) and 



fd^xf{x)S{x,K) 
J d'^xS{x, K) 



(8) 



(9) 



an average with the emission function. 



2.6 Fit results — mj-spectra and HBT radii 

The model is fitted to the recent SPS data by NA49 [18], hence there are small differences to the results 
of [20] used in [4-6, 8, 9] As an example, we show the resulting transverse mass spectra for tt^ and 
in Fig. 12 (left) and the rapidity distribution for n~ (right) compared with the NA49 data. 

There is a contribution to the 7r~ transverse mass spectrum coming from the vacuum decay of resonances 
after r/ . The model is able to calculate the magnitude of this contribution using statistical hadronization 
but not its nit distribution, therefore it is left out when we compare with data. This explains the 
mismatch at low transverse mass and that the integrated spectrum does not yield the multiplicity 
indicated by the data. The same is visible in the spectrum, albeit less pronounced. Apart from this 
aspect, the model describes the transverse mass spectra well. 

The fact that particles inside the fireball volume are characterized by a thermal momentum distribution 
leads to a smearing of the box-profile in rapidity when particle emission is calculated, therefore the 
resulting 7r~ distribution in rapidity docs not exhibit a sharp cutoff. The deviation from the data at 
forward rapidities is caused by the simplifying assumption of a homogeneous rapidity density, it is not 
a failure of the framework as such and indicates that this approximation should not be used for y > 1.5. 
We illustrate this by replacing the box density distribtion in rapidity with a Woods-Saxon distribution 
with a skin thickness parameter Ay = 0.37 which dramatically improves the agreement with the data. 
We have checked that this choice does not significantly alter the transverse mass spectra at y = 0. 

The resulting HBT correlation parameters are shown in Fig. |21 In general, the agreement with the data 
is good, however the low kt region seems to be underestimated. Note however that our calculation does 
not contain the model independent correction terms which tend to increase the correlation radii for small 
momenta (see our results in [10]) and that the experimental data do not include systematic errors. 
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Figure 2: Left panel: Transverse mass spectra of tt^ and measured by NA49 [18] compared with the 
model results. Right panel: The rapidity distribution of tt^ measured by NA49 [18] (defining as the 
rapidity of the center of mass system) compared with the model result (solid) and a calculation replacing 
the box density profile by a Woods-Saxon distribution with skin thickness Ay = 0.37 (dashed). 



2.7 The fireball evolution 

The model for 5% central 158 AGeV Pb-Pb collisions at SPS is characterized by the following scales: 
Initial longitudinal velocity Wq ™* — 0.54c, equilibration time tq = 1 fm/c, initial temperature Tq — 300 
MeV, duration of the QGP phase tqgp = 6.5 fm/c, duration of the hadronic phase Thad = 8.5 fm/c, 
total lifetime Tf — tq = 15 fm/c, r.m.s radius at freeze-out i?""^ = 8.55 fm, transverse expansion velocity 
v±f — 0.57c and freeze-out temperature Tp — 100 MeV . Fig.Qlshows the resulting temperature evolution 
in proper time and compares it with the previous fit results based on older data used in [4-6,8,9]. 

The difference between the two fits is not large — the new data require a slightly larger Uq'°"' , leading to 
a reduced initial temperature, but this change is well within the model uncertainties. A more pronounced 
difference shows up in the late hadronic evolution where a larger is required, leading to faster cooling 
and a reduced lifetime in the new fit. Surprisingly, this amounts to an increase of the four- volume of 
the hadronic phase of about 5%, the shorter time duration is slightly outweighed by a faster volume 
increase. 

Based only on a fit to hadronic observables, we thus find a system which is characterized by a sizeable 
initial compression of matter with subsequent re-expansion (wg'™' < u^™"*). This has several important 
consequences: First, any estimate of the initial temperature or total lifetime based on a boost-invariant 
expansion such as the well-known relation i?iong = TfiTp/'n^tY^'^ connecting i?iong with the breakup 
time [21] does not apply. Instead, a system undergoing accelerated longitudinal expansion is initially 
more strongly compressed (implying a larger initial energy density and temperature) and takes more time 
to expand to a given volume while arriving at the observed longitudinal flow. Hence, the four-volume of 
such a system is also larger than in the standard boost-invariant expansion scenario. We argue later that 
the difference between the two cases is experimentally accessible by measuring electromagnetic probes. 

For the discussion of some observables, we require the fireball evolution for other than 5% central 
collisions. In this case, we scale the total entropy with the number of collision participants and reduce 
the initial radius of the system to reproduce the smaller overlap area of the nuclei while neglecting the 
angular asymmetry. Keeping the parameter c± at its value, the reduced entropy leads to an earlier 
freeze-out and hence naturally reduces transverse flow. For the initial longitudinal motion, we linearly 
interpolate between the value found for central collisions and the observed rapidity loss in p-p collisions 
(where no re-expansion occurs) in the number of collision participants and refit Cz accordingly. 
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Figure 3: HBT Correlation radii obtained by the CERES collaboration [19] as a function of trans- 
verse pair momentum kt compared with the model results. No systematic errors are included in the 
experimental errorbars. 



3 Statistical hadronization at the phase boundary 

In order to connect the fireball entropy with the measured particle spectra, we have to specify how 

a hadronic system at the phase boimdary is created. Statistical models arc extremely successful in 
describing the measured ratios of different hadron species for a range of collision energies from SIS to 
RHIC (see e.g. [22-25]). 

The basic assumption of statistical hadronization is that the particle content of the fireball can be found 

by considering a system of (non-interacting) hadrons in chemical equilibrium at Tc, described by the 
grand canonical ensemble, which subsequently undergoes decay processes. This ensemble is characterized 
by the temperature Tc, the baryochemical potential fiB and the strange potential /Ug. In the present 
model framework, we can calculate all these parameters from the fireball evolution and obtain hadron 

ratios parameter-free. 

Employing the grand canonical ensemble, the density for each particle species rii is given by 

2772 J, exp{ [E^ip) - ^l,]/Tc} ± 1 ■ ^ ^ 

Here, di denotes the degeneracy factor of particle species i (spin, isospin, particle / antiparticle) , the 
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Figure 4: Temperature evolution of the fireball for SPS 158 AGeV 5% centeral Pb-Pb collisions, shown 
are the old fit used in [4-6,8,9] and the scenario based on recent data presented in this work. 



+ (-) sign is used for fermions (bosons) and Ei{p) — \/rn^ + p^. rrii stands for the particle's vacuum 
mass. We use a value of 170 MeV for the critical temperature Tc as determined in lattice calculations 
for two light and one heavy flavours [2,3]. The chemical potential fii = fJ-sBi — fisSi takes care of 
conserved baryon number Bi and strangeness Si for each species and we neglect a (small) contribution 
—^ii^lf coming from the isospin asymmetry. The conservation of the number of participant baryons 
A^part and the requirement of vanishing net strangeness uniquely determines and /ig for a given 
volume V = V{Tc)- This volume can be calculated as V{Tc) = Stot/s{Tc) with s{Tc) determined by 
the EoS in the partonic phase, thus linking the model calculation again with lattice results. 

We include all mesons and mesonic resonances up to masses of 1.5 GeV and all baryons and baryonic 
resonances up to masses of 2 GeV in the model and calculate their decay into particles which are stable 
or long-lived as compared to the fireball, such as tt, K, 77, N, A, E and SI. In order to account for repulsive 
interactions between particles at small distances we assume a hard core radius Rc of 0.3 fm (based 
on results for p-p collisions [26]) for all particles and resonances and correct for the excluded volume 
self-consistently. Particle and resonance properties are taken from [27]. For resonances with large width, 
we integrate Eq. (|10|l over the mass range of the resonance using a Breit-Wigner distribution. 

10.00 r , 
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Figure 5: Hadron ratios in the statistical hadronization model (dashed bands) as compared to experi- 
mental results (filled circles) for SPS, 158 AGeV central Pb-Pb collisions [28-38]. 
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The result (Fig. O is in good agreement with the experimental data [28-38] . The ratios are (for given 
Tc) determined by the volume which in turn determines ^b- We may infer from this agreement that the 
value of Tc as measured on the lattice is, in combination with the statistical hadronization approach, 
compatible with the data and that the EoS, via s(Tc), leads to a reasonable volume at the phase 
transition. Therefore, hadron ratios are insensitive to the spacetime expansion parametrized in the 
model but test mainly properties of the EoS. 

The main impact of the results on the fireball evolution is via the EoS: Decay processes of high- lying 
resonances created at Tc lead to on overpopulation of pion phase space during the evolution in the 
hadronic phase when strong decay processes are still in equilibrium with their back reaction but weak 
decays are not. This effect has a crucial influence on the EoS and we incorporate it by allowing for a 
finite pion chemical potential iJ-Trir) that rises to good approximation linearly in proper time from at 
r = rc to 110 MeV at r = (cf. [39]). 



4 Electromagnetic observables 
4.1 Dilepton emission 

The lepton pair emission rate from a hot domain populated by particles in thermal equilibrium at 
temperature T is proportional to the imaginary part of the spin-averaged photon self-energy, with these 
particles as intermediate states. The thermally excited particles annihilate to yield a time-like virtual 
photon with four-momentum q which decays subsequently into a lepton-antilepton pair. The differential 
pair production rate is given by [40] 

' in.n(,.T)^^J|^, (11) 



where a — e^/An, (3 = 1/T, the four- velocity of the emitting volume element and we have neglected 
the lepton masses. We have defined Ii^{q) — —11^^/3 and introduce the averaged photon spectral function 
R{q,T) — (1277/(7^) lm1J{q,T). Here 11^^ denotes the trace over the thermal photon self-energy which is 
equivalent to the thermal current-current correlation function 

n^,(g,T) = tjd^x e'i^Tj^{x)j,{()))p, (12) 

where is the electromagnetic current. Ea. (|ll() is valid to leading order in the electromagnetic interac- 
tion and to all orders in the strong interaction. The information about the strong interaction dynamics 
is now encoded in the spectral function. For its computation, we have to use different techniques for 
dealing with partonic and hadronic degrees of freedom. 



4.1.1 Spectral function in the QGP phase 

As long as the thermodynamically active degrees of freedom are quarks and gluons, the timelike photon 
couples to the continuum of thermally excited qq states and subsequently converts into a charged lepton 
pair. The calculation of the photon spectral function at the one-loop level is performed using standard 
thermal field theory methods. The well-known leading-order result for bare quarks and gluons as degrees 
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of freedom is [40]: 

f=u,cLs \ ^ / 

X 1^1 + 2 

where q = (g^jq) is the four-momentum of the virtual photon, e/ the quark electric charge, ^{q^) = 

yl ^ and ruf the quark mass of flavour /. This result, however, is modified by perturbative 

corrections in as that take into account interactions in the plasma. 

We compute these within our quasiparticle description [12, 13]. Here, the degrees of freedom in the par- 
tonic phase acquire an effective temperature dependent mass m(T). Furthermore, the effect of confine- 
ment near the phase transition is parametrized by a reduction factor C(T) multiplying the distribution 
functions (see references for details). 

It is straightforward to insert thermal masses instead of bare quark masses into Eq. 1131 C(T) is equivalent 
to a fugacity factor, it describes the reduced occupation of available states for quarks and antiquarks 
in the medium due to confinement. Since Eg. 1131 describes the annihilation of a qq pair into a virtual 
photon, the probability of finding a quark and antiquark in the medium is reduced by a factor C(T)^ in 
this approximation. 




4.1.2 Spectral function in the hadronic phase 

Below Tc, confinement sets in and the effective degrees of freedom change to colour singlet, bound qq 
or qqq (qqq) states. The photon couples now to the lowest-lying dipole excitations of the vacuum, the 
hadronic = 1~ states: the p, cu and (p mesons and multi-pion states carrying these same quantum 
numbers. There is a considerable uncertainty in the calculation of properties of hadronic matter near 
the phase boundary. In order to minimize the theoretical uncertainty, we compare several different 
approaches: 

The first approach (referred to as 'chiral model' in the following) is based on [41,42,44]: The elec- 
tromagnetic current-current correlation function can be computed from an effective Lagrangian which 
approximates the SU (3) flavour sector of QCD at low energies. An appropriate model for our purposes 
is the Improved Vector Meson Dominance model combined with chiral dynamics of pions and kaons as 
described in [43]. Within this model, the following relation between the imaginary part of the irreducible 
photon self-energy Imll and the vector meson self-energies Ilv{q) in vacuum is derived: 

I-n(g) = >^ \Fv{q)\ , Fyiq) = ^2 _ + ,i^nv(g) ' ^'^^ 



where my are the (renormalized) vector meson masses, gy is the jV coupling and g is the vector meson 
coupling to the pseudoscalar Goldstone bosons tt*, 7r° and K^, Ea. l|14|l is valid for a virtual photon 
with vanishing three-momentum q. For flnite three- momenta there exist two scalar functions 11^ and 
ILt, because the existence of a preferred frame of reference (the heat bath) breaks Lorentz invariance, 
and one has to properly average over them. However, for simplicity we approximate the problem by 
taking the limit |q| — )■ and test this approximation later in a different approach. 

Finite temperature modifications of the vector meson self-energies appearing in ea. (|14|l are calculated 
using thermal Feynman rules. The explicit calculations for the p- and (/)-meson can be found in [41]. 
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The thermal spectral function of the w- meson is discussed in detail in [42]. The p mass is not shifted 
explicitly in the approach. 

For the evaluation of finite baryon density effects which are relevant at SPS conditions, we use the results 
discussed in [44]. There it was shown that in the linear density approximation, Xly is related to the 
vector meson - nucleon scattering amplitude. In the following, we assume that the temperature- and 
density-dependences of Hv factorize. This amounts to neglecting contributions from matrix elements 
such as (7rA^|Tj^(a;)j'^(0)|7r7V) describing nucleon-pion scatterings where the pion comes from the heat 
bath. 

Our second approach [9] (referred to as 'mean field model' in the following) is based on a different idea: 
The method of thermofield dynamics [45] is used here to calculate the state with minimum thermody- 
namic potential, at finite temperature and density within the mean field sigma model with a quartic 
scalar self interaction. The temperature and density dependent baryon and sigma masses are calculated 
self-consistently. The medium modification to the masses of the lo- and p-mesons in hot nuclear matter 
including the quantum correction effects are then calculated in the relativistic random phase approxi- 
mation. The decay widths for the mesons are calculated from the imaginary part of the self energy using 
the Cutkosky rule. 

In a third approach (referred to as 'tt — p-model' in the following) [46] , we solve truncated Schwinger- 
Dyson equations of the tt — p system in a self-consistent resummation scheme, neglecting the effects of 
other hadrons. This approach takes into account the finite in-medium damping width of the pion which 
contributes substantially to the broadening of the p-meson. As an additional benefit, we can use this 
model to study the effects of neglecting the q dependence of the spectral function in chiral model. 



4.1.3 The measured spectrum 

The differential rate of Ea. Hll() is integrated over the space-time history of the collision to compare 
the calculated dilepton rates with the CERES/NA45 data [47] taken in Pb-Au collisions at 158 AGeV 
(corresponding to a cm. energy of ^/snn ~ 17 AGeV). The CERES experiment is a fixed-target 
experiment. In the lab frame, the CERES detector covers the limited rapidity interval rj = 2.1 — 2.65, 
i.e. A?7 = 0.55. We integrate the calculated rates over the transverse momentum pT and average over 
rj. The resulting formula for the space-time- and p-integrated dilepton rates is 

Tf R{t) >)/™„t(T) OO 

d^N 2nM [ f f sinh(C - 1) , f, diV(r(T), Af, ry(r?,),pT) , . 

(15) 

where the function Acc(M, 77, pt) accounts for the experimental acceptance cuts specific to the detector. 
In the CERES experiment, each electron/positron track is required to have a transverse momentum 
Pt > 0.2 GeV, to fall into the rapidity interval 2.1 < 77 < 2.65 in the lab frame and to have a pair 
opening angle Oee > 35 mrad. Finally, for comparison with the CERES data, the resulting rate is 
divided by dNch/dr], the rapidity density of charged particles. 

In addition to the thermal emission of dileptons, we also consider dileptons from decays of vector mesons 
p,uj,ip after the kinetic decoupling of the fireball. We calculate their yields using the statistical hadroniza- 
tion model outlined in sectionOand use their vacuum spectral functions to calculate the resulting dilepton 
yield. 
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4.1.4 Results 



The invariant mass spectrum of dileptons is a Lorentz-invariant quantity and hence independent of flow 
and roughly proportional to the radiating 4-volume in the hadronic phase. This explains the negligible 
differences between the old and the new analysis in spite of the fact that the refitted evolution scenario 
exhibits stronger transverse flow and has a reduced lifetime. 

The calculated spectrum is in fair agreement with the data for all three spectral functions (see Fig.O; 
the mean field model seems superior to the chiral model in the description of the invariant mass re- 
gion below 300 MeV, however, the general magnitude of dilepton emission is well reproduced by both 
spectral functions. As Fig. (right) indicates, a purely pionic medium is also able to account for a 
sufficient broadening of the p meson if the pion damping width is taken into account in a self-consistent 
resummation scheme instead of using a perturbative expansion. 

Neglecting the momentum dependence in the tt — p model to overestimates the data in the low invariant 
mass region. This is the same trend seen in the chiral model and suggests that including the momentum 
dependence there would probably also improve the agreement with the data. 




0.5 1 1.5 

Invariant mass [GeV] 



• CERES data 158 AGeV 
— Self-consistent 7l-p model 
■ ■ neglecting 3-momentum dep. 




(1.5 1 1.5 

Invariant mass [GeV] 



Figure 6: Left panel: The dilepton spectrum compared with SPS 158 Pb-Au results [47] (dots) using 
both the Chiral model (dashed) and the Mean Field model (solid) spectral function. Right panel: The 
dilepton invariant mass spectrum in the self-consistent tt — p model (solid) and in a calculation neglecting 
the 3-monientuni dependence (dotted). Note that although the Mean Field model spectral function does 
not contain the (j) meson and the tt — p model contains neither lo nor 0, vacuum decays of these mesons 
are nevertheless included and lead to distinct peaks in the spectrum. In both plots we did not fold our 
result with the finite energy resolution of the detector, this explains the apparent disagreement of the 
sharp (j) and lo vacuum decay peaks with the data. 

To gain additional insight into the underlying physics, we convolute the chiral model result with the 
finite energy resolution of the detector and indicate the individual contributions to the measured rate 
separately. This is shown in Fig. |7| The decomposition for the Mean Field model is very similar. 

The dominant contribution to the complete rate in the invariant mass region below 1 GeV comes from 
thermalized hadronic matter, and here by far the strongest contribution from the p. For both the Chiral 
and the Mean Field model calculation, the inclusion of finite baryon density effects (which manifest 
themselves in N-p scattering processes) is crucial for achieving agreement with data [4,9]. Both calcula- 
tions take these into account. The success of the n — p model suggests, however, that even for vanishing 
baryon density the p-meson may be sufficiently broadened to explain the data. 
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Figure 7: Left panel: Dilepton emission rate in the Chiral model folded with the energy resolution of 
the detector (thick soUd) compared with SPS 158 AGeV Pb-Au results (solid circles). Indicated are 
the thermal contribution from the QGP (thin solid), contributions from initial hard Drell-Yan processes 
(dash-dot) and vacuum decays of mesons, the so-called cocktail (dashed). Right panel: The total thermal 
emission rate from hadronic degrees of freedom (solid), decomposed into contributions from p (dashed), 
uj (long dashed) and (j) (dotted). 



The QGP contribution is only visible above 1 GeV invariant mass, however here the large errors do 
not allow any firm conclusion about its magnitude. As is apparent from Fig. \7\ the QGP contribution 
is unable to explain the enhancement of the dilepton spectrum below 600 MeV invariant mass. Thus 
a dilepton measurement at SPS is dominantly sensitive to the in-medium mass modifications of vector 
mesons. While those are interesting problems in their own right, it implies that the dilepton measurement 
at CERES does not directly help to clarify the question if the quark-gluon plasma has been formed. 

The importance of thermal contributions from the hadronic phase serves as a reminder that interactions 
do not cease to be important after the phase transition. Any type of model assuming kinetic freeze-out 
at the phase transition would have to find a different mechanism of producing the required amount of 
dilepton radiation below the p peak — clearly, the QGP contribution alone is insufhcient to explain the 
enhancement. On the other hand, the fact that the hadronic spectral functions we have investigated 
lie on the upper and lower bound of the error bars below 0.5 GeV invariant mass indicates that the 
hadronic contribution (dictated by the four-volume of the radiating matter) is of the correct order of 
magnitude. If the radiating 4- volume were less than half of the value in the present model, the calculated 
hadronic emission spectrum would not be sufficient to explain the data. This strongly disfavours large 
freeze-out temperatures and a short lifetime of the hadronic phase (although for sufficient transverse 
flow, a large emitting 4- volume could still be achieved in spite of a small lifetime; the measured hadron 
spectra, however, rule out such a combination). 



4.2 Photon emission 
4.2.1 The emission rate 

For the emission rate of direct photons, we use the complete 0{as) calculation [48-54] in the form of 
the parameterization provided in [54]. Using arguments analoguous to the ones in the case of dilepton 
emission, we incorporate the effect of the quasiparticle description into this equation by a factor C{T)^. 
For the photon emission rate from hadronic matter, we use a parametrization of the rate from a hot 
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hadronic gas taken from [55]. 



4.2.2 The emission spectrum 

The spectrum of emitted photons can be found by folding the rate [54] with the fireball evolution. In 
order to account for flow, the energy of a photon emitted with momentum fc^ = (fc(,kt,0) has to be 
evaluated in the local rest frame of matter, giving rise to a product k^Ufj, with u^{ris,r, t) the local flow 
profile as decribed in section 13 For comparison with the measured photon spectrum [56], we present the 
differential emission spectrum into the midrapidity slice y = 0. The resulting photon spectrum is shown 
in Fig. IHl (left). 




k| [GeV] k| [GeV] 



Figure 8: Left panel: Thermal photon spectrum for 10% most central Pb-Pb collisions at SPS, 158 
AGeV Pb-Pb collisions, shown are calculated rate (total, contribution from QGP and hadronic gas) and 
experimental data [56]. Right panel: The calculated photon spectrum (solid), the spectrum assuming 
boost-invariant longitudinal expansion (dashed) and assuming formation times 0.5 fm/c (dash-dotted) 
and 3 fm/c (dotted) in comparison with the data [56]. 

The overall agreement with the data is remarkably good. We observe that the relatively large transverse 
flow in the hadronic phase in the present scenario leads to a rather flat contribution with a shape 
similar to the QGP contribution but smaller. For the standard choice tq — 1 fm/c the calculation 
seems to underestimate the data slightly in the high momentum tail, but yields a good description over 
the whole range 1.5 GeV < fct < 3 GeV. The errors allow for a prompt photon contribution of about 
half of the magnitude of the thermal yield for the whole momentum range. Note that the spectrum is 
almost completely saturated by the QGP contribution — the hadronic contribution is about an order of 
magnitude down. 

In order to study the influence of the equilibration time tq and hence the initial temperature, we present 
model calculations for the choices tq = 0.5 fm/c and tq = 3 fm/c in fig. |S1 (right). Clearly, a late 
equilibration fails to account for the data in the momentum region between 2 and 2.5 fm/c where prompt 
photons are not expected to play an important role (cf. the discussion in [6]). On the other hand, an 
equilibration as early as 0.5 fm/c (corresponding to an initial temperature of 370 MeV) produces too 
many photons. The direct photon spectrum thus constrains the initial equilibration time at SPS energies 
to be close to 1 fm/c. 

We may finally use the photon emission rate to test the validity of our fireball model. Tentatively 
assuming a standard boost-invariant expansion pattern with ?7q'°"* = 77^''°"', we calculate the resulting 
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photon emission. This is also shown in Fig.|Sl right panel. The rapid cooling and the consequently short- 
lived QGP phase for such a scenario lead to a severe disagreement with data that cannot be compensated 
for by any choice of equilibration time above 0.1 fm/c. Any initial state photon emission is bound to 
fill the larger kt region of the emission spectrum and cannot properly compensate for the reduction of 
radiating four-volume at moderate temperatures (below 250 MeV) in the boost-invariant scenario. We 
conclude that longitudinal compression and re-expansion of matter with all its consequences is essential 
for the calculated thermal photon yield to be compatible with the data and that the analysis indicates 
a rather quick equilibration. 

5 Charmonium suppression 

It is often argued that the breakup of bound cc states immersed into hot and dense matter is a good signal 
for the creation of a QGP. We investigate this question within the framework of our fireball evolution 
model. We refer to J/^ and the excited states ^' and Xc generically as ^ in the following. Details 
of our treatment of the excited states can be found in [8]. Our main assumption is that charmonia 
themselves are not thermalized. This is motived by the observation that the mass scale of charm quarks 
is well above the temperature scale even for the initial conditions. The thermal production of cc pairs 
is therefore negligible and the momentum transfer from produced charm quarks to light constituents of 
the heat bath is not very efficient. 

Based on this assumption, we will treat the formation and breakup of bound charmonia as an off- 
equilibrium process and use kinetic theory to calculate the evolution of the ^' density. This corresponds 
to the treatment in [57]. We will here limit the discussion to ^E* dissociation in the QGP phase in the 
following because the density of scatterers in the medium as a function of temperature in the hadronic 
phase is two orders of magnitude smaller than in the QGP phase [8]. This is illustrated in Fig. |^ which 
shows the particle density. (Fig. ^ 
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Figure 9: Density of medium constituents as calculated in the quasiparticle model [12, 13] as a function 
of temperature. 



5.1 Charmonium production in a nuclear environment 

We start by parametrizing charmonium production in p-p collisions. For later use we need the pt 
spectrum of ^E* at mid-rapidity. In the following we assume a Gaussian form for the p^-dependent 
part, with width parameter A = 1 GeV/c. The rapidity modulation can be inferred from the relation 
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da/dy ^ xig(xi) X2g{x2) where xg{x) ~ (1 — x)^ is the gluon distribution in the proton and xi^2 = 
(mif/y^) exp(±y). For the overall normalization we use the paramctrization for the total charnionium 
production cross section [58] 

a*(s) = 2ao(l-m*/Vi)", (16) 

where ctq — 1-28 and n — 12. 

We now consider nuclear effects, starting with the simpler case of proton-nucleus {pA) collisions. It has 
been shown that the experimental results on charmonium (^) production can be explained using [58] 

o^A^^j d^^TA{h)ST{h) (17) 

for the total production cross section. The factor 



1 — exp 



Sa (b) = -77-, (18) 

^f^TAib) 

is the survival probability for "if to escape the nucleus without being dissociated. It includes the effective 
absorption cross section crj^'j^, a quantity of the order of 3 mb for mid-rapidity ^fs as measured at Eiab = 
800 GeV at Fermilab, while it amounts to 5 — 7 mb for mid-rapidity 4's as measured at Eiab = 158 — 200 
GeV at the SPS. The absorption cross section parametrizes various poorly known effects, with varying 
importance depending on the collision energy. Among these effects are the presence of color degrees of 
freedom in the dynamics of colliding nucleons, initial state parton energy loss and coherence length and 
shadowing effects. A common property of all of the above is the linear dependence on the path length, 
at least to leading order. Using eq. H18() can therefore be justified, provided a suitable re-scaling and 
re- interpretation of cr^^ ct^n and Sa'^% ^a^^ is done. 

When considering 5" production in nucleus-nucleus (AB) collisions, one can estimate the cross section 
for a given impact parameter by generalizing eq. ifTTjl . Neglecting effects of the medium produced in the 
collisions (which will be discussed in the following) one obtains 

where nuclear effects are included in the suppression function 

STB^ib) = T^hib) I d's Ta{s) Sr^'is) TsiS) 5r^(5). (20) 
We choose tr^^ = 5 mb at the SPS energy y/s = 17.3 GeV conforming with the pA measurement 



5.2 Kinetic description of \1/ breakup 

We describe the suppression of charmonia by using kinetic theory. Neglecting the possibility of charmo- 
nium formation in the medium for SPS energies, we can derive the expression [8] 

^ = -E((->-))^»^-A.(r)^ (21) 

for the evolution of the charmonium rapidity density. Here, crj) denotes the dissociation cross section for 
collisions with medium constituent n, pn its density and the sum runs over all possible constituents. The 
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double brackets indicate an average over the momenta of the initial-state particles ^ and n except the 
^ rapidity, we want to compute the final '5 rapidity distribution at mid-rapidity, not the whole yield. 

The ingredients required to solve this equation are the initial distribution of charmonium discussed in 
the previous section, the evolution of the medium and the functional form of gd- 

We use the result of Bhanot and Peskin [59,60] 



for the gluon dissociation process of a heavy quarkonium. It is a function of the gluon energy lo in 
the rest frame of the quarkonium and contains the threshold energy e and the mass scale /i, related to 
the heavy quark mass. The threshold energy is related to the binding energy cq by the condition that 
s = {p + fc)^ > 4m^ which implies that e = + eQ/(2rniif). In the following we assume eg = 780 MeV 
for the binding energy and /Xc = 1.95 GeV for the mass as in [60] to fit the mass values of the first two 
levels ( J/'0 and tp') of the charmonium system. 

5.3 Results 

We combine the building blocks described above to compute the suppression as a function of impact 
parameter h. The result is divided by the number of binary collisions A'coii, which provides the centrality 
dependence of the Drell-Yan cross section. The normalization is fixed at A^coii = 53.5 to match, at 
large 6, the experimental value in p-p collisions and the ratio / DY as function of h is converted into 
a function of the measured transverse energy assuming ET{b) — iVpart(&)- The quantity et = 0.297, 
which represents the amount of produced transverse energy per participant, is used in order to describe 
correctly the total inelastic (minibum-bias) cross section as function of centrality. In this way we arrive 
at the results plotted in Fig. IIUI 

The agreement with data is quite remarkable but deserves some comments. First, the curves end at 
Et « 110 GeV, which corresponds to 6 = 0. To go beyond this value it is necessary to include eff'ects 
of fluctuations, which are in principle straightforward to address [61]. However, as the analysis in [62] 
suggests, fluctuations do not alter the results below this value substantially, therefore such an extension 
of the model would not yield new information about the scales of the spacetime expansion of the fireball, 
the main issue under investigation in this paper. 

The agreement of our curve with the data down to Et = 40 GeV is remarkable given the fact that we 
use rather simple scaling laws [8] to model the impact parameter dependence of the fireball evolution. 
On the other hand, the deviation from the date for Et < 40 GeV is not surprising, because the thermal 
equilibrium assumption is bound to fail for very peripheral collisions. 

It is reassuring that the model is able to describe the data over a wide range, however there are still con- 
siderable uncertainties regarding the excited ^' states, the magnitude of the 'ordinary' nuclear suppression 
and the magnitude of the production cross section for nucleus- nucleus collisions (cf. [8]). Therefore, we 
refrain from deriving constraints for the model framework from this observable and note only that the 
same early evolution which describes the photon data also agrees well with the charmonium suppression 
data. 



6 Summary 

Under the assumption that a thermalized system is created in 158 AGeV Pb-Pb collisions at SPS, we 
have analyzed a schematic model based on the thermodynamic response of matter to a volume expansion. 




(22) 
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Figure 10: Result at SPS energy for the ratio ^/DY as function of the transverse energy. The dashed 
curve includes only nuclear effects, while the full line is the complete result including gluon dissociation 



This model provides the link between the EoS as obtained in lattice QCD and measured observables. We 
have derived essential evolution scales by the requirement that the model should reproduce the measured 

properties of hadron emission. 

Various other sets of observables provide unique opportunities to verify these essential scales. Discussing 
statistical hadronization, we find that the crucial quantity in our approach is given by s{Tc), the entropy 
density at the phase transition. This quantity is an input from lattice QCD and agrees nicely with the 

observed hadon ratios. 

A calculation of dilepton emission shows that most of the observed dilepton excess above the 'cocktail' 
of final state meson decays can be explained by radiation from a thermalized hadron gas. The main 
contribution here stems from a strongly broadened p-meson. Thus, we find evidence for ongoing inter- 
actions even below the phase transition temperature. The total amount of radiating space-time-volume 
in our model description is in good agreement with the measured amount of radiation. 

The measurement of high-momentum photons gives complementary information: The large momentum 
scale makes this observable sensitive to the initial state of the collision. Again, we find good agreement 
with the data and explicitly demonstrate that a standard boost-invariant expansion does not yield 
a satisfactory description. Therefore, initial longitudinal compression and subsequent expansion, an 
essential point of our framework, seems to be realized at SPS conditions. 

This observation is further strengthened by a calculation of charmonium suppression, which we also find 
to be predominantly sensitive to the initial state. Although uncertainties in the calculation are not small 
and many simplifications have been made, the agreement with data is certainly reassuring. 

In summary, we have presented a thermal description of relativistic heavy-ion collisions at the CERN 
SPS for 158 AGcV Pb-Pb which is highly consistent and leads to a good description of a large set of 
different observable quantities. While this docs not provide a conclusive proof for the creation of a 
thermalized system (and hence a QGP), it constitutes certainly strong evidence for it. It remains to 
be investigated if the results can be reproduced in a thermal microscopical transport description (like 
hydrodynamics) and if a non-thermal framework is able to describe the experimental results as well or 
if thermalization is required by the data. 
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